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Dh: abstract 
o 

^ \ Dwarf irregular galaxies are unique laboratories for studying the interaction 

between stars and the interstellar medium in low mass environments. We present 
the highest spatial resolution observations to date of the neutral hydrogen content 
^ ' of the Local Group dwarf irregular galaxy WLM. We find that WLM's neutral 

d \ hydrogen distribution is typical for a galaxy of its type and size and derive an 

HI mass of 6.3 ± 0.3 x 10^ Mtn for WLM. In addition, we derive an HI extent for 



WLM of ~ 30', which is much less than the 45' extent found by lHuchtmeier et al. 



( 1l98ll ). We show that the broken ring of high column density neutral hydrogen 
surrounding the center of WLM is likely the result of star formation propagating 
out from the center of the galaxy. The young stars and Ha emission in this galaxy 
are mostly correlated with the high column density neutral hydrogen. The gap 
in the central ring is the result of star formation in that region using up, blowing 
out, or ionizing all of the neutral hydrogen. Like many late-type galaxies, WLM's 
velocity field is asymmetric with the approaching (northern) side appearing to 
be warped and a steeper velocity gradient for the approaching side than for the 
receding side in the inner region of the galaxy. We derive a dynamical mass for 

0- 



WLM of 2.16 X 10^ M^. 
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Subject headings: galaxies: dwarf, galaxies: irregular, galaxies: ISM, galaxies: 
kinematics and dynamics. Local Group, galaxies: individual (WLM) 



Introduction 



The interstellar medium of dwarf irregular galaxies provides important information 
about how stars and gas interact in low mass environments. Feedback from massive stars 
is the dominant process that affects the interstellar medium in these galaxies. Their low 
mass makes it easier for violent episodes of star for mation to disrupt the gala xy's gas, pos- 
sibly even expelling it from the galaxy's potential (IMac Low &: Ferraralll999l ). The lack of 
differential rotation in dwarf irregulars contributes to longer lifetimes for features like holes 
in the interstellar medium created by star formation. These longer hole lifetimes allow one 
to probe the interaction of massive stars and the interstellar medium over longer timescales 
than possible in other types of galaxies. 

As part of an ongoi ng survey of dwarf irregular galaxies with appare ntly extended neu- 
tral hydrogen emission (jWilcots &: Hunterl l2002l : iHunter fc Wilcotsll2002l ). we have obtained 
observations of the neutral hydrogen, ionized hydrogen, and stellar content of an isolated 
Local Group dwarf irregular galaxy called WLM. Its nearest neighbor, the Cetus dwarf 
sphe roidal galaxy, is appr oximately 200 kpc away, assuming a distance to Cetus dSph of 775 
kpc (jWhiting et al.lll999l ). Half of WLM's stars, by mass, were formed over 9 Gyr ago, and 
since then it has been forming stars at a leisurely rate of 1 — 2 x 10"'^ M0 yr~^ (IDolphin 
2OOOI ) . In the last Gyr or so, there has been a spurt of star formati o n, with new stars being 
forme d at three times the average rate in the past (jDolphirul2000l ). lYoungblood fc Hunter 
(119991 ) identified 30 Ha regions in WLM, which is consistent with its recent spurt of star 
formation. Like many dwarf i rregulars, WLM has a low metallicity; its oxyg en abundance 
is 12 + logO/H = 7.78-8.23 Jskillman et al.l Il989l : iHunter fc HoffmanI Il999h . where 8.7 is 



solar. Given its metallicity, it is not surprising that iTaylor fc KleinI (I2OOII ) did not detect 
any CO emission in this galaxy, since the CO to H2 ratio is thought to decrease significantly 
for low metallicity ga laxies. WLM also has a globular cluster with an age of 14.8 ± 0.6 Gyr 
on its western edg e (IHodge et al.lll999l ). which is very unusual for a faint dwarf irregular 
teillett et all 12002! ). 



We present the highest spatial 
of WLM. Previous studies include 
locity resolution: 1.65 km s~^ V 



resolution study to date of the neutral hydrogen content 



Jackson et al.l (120041) f spatial resolution: 39" by 161"; ve- 
Barnes &: de Blo < (12004 ) (spatial resolution: 15.5'; velocity 



resolution: 1.12 km s ^ ), and iHuchtmeier et al.l (Il98ll ) (spatial resolution: 9'; velocity reso- 
lution: 11 km s~^ ). WLM's declination makes it difficult to observe; it is low enough in the 
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sky for Very Large Array's (VLA) beam to be elongated and is too near the celestial equator 
to get good uv coverage with the Australia Telescope Compact Array (ATCA). We have ob- 
tained vlaQ observations of WLM in the B, C, and D configurations. Although our beam is 
elongated in the north-south direction due to the low declination, the high spatial resolution 
of these observations allows us to investigate in detail the HI content of this galaxy. We 
have complemented these observations with UBV and Ha imaging to gain a more complete 
picture of the interaction between the stars and the interstellar medium in WLM. Our data 
are discussed further in Section [2J In Section [31 we use our data to examine the properties 
of the total HI content of WLM. We examine the connection between the small-scale struc- 
ture of WLM's interstellar medium and its stars in Section |H Section [5] is a detailed look 
at the rotation curve of WLM (associated mass modeling is shown in Appendix Our 
conclusions are summarized in Section O 



2. Data 
2.1. VLA observations 



Previous observations of WLM (IHuchtmeier et al.lll98ll ) indicated that the HI extent of 



WLM is larger than the primary beam of the VLA at 20 cm (45' versus 30'), so we mapped 
WLM in D array as a 3 by 3 mosaic with pointing centers 15' apart to make sure that we 
would detect the total extent of the large scale structure. The resulting square field of view 
is 60' (at the half-power points). These observations were taken on 2001 October 19. The 
data were on-line hanning smoothed to give a final channel width of 2.57 km s~^ . The flux 
calibrator (0137+331) was observed at the beginning and end of the observations. During 
the observations, we cycled through the mosaic positions, observing each for 100 s per cycle. 
We observed the phase calibrator (2321-163) once every two cycles through the mosaic. The 
total integration time for each pointing was 28 minutes. The effective total integration time 
for positions near the center of the mosaic was between 128 and 140 minutes due to the 
overlap between the pointings. 

To obtain information on the small scale HI structure of WLM, we observed it using 
the VLA in B and C configurations. The B array observations were taken on 2002 July 9 
and 12 and consisted of a single pointing centered on WLM's position. The fiux calibrator 
(0137-1-331) was observed for 15 minutes at the end of each day and the phase calibrator 



^The VLA is part of the National Radio Astronomy Observatory, which is a facihty of the National 
Science Foundation operated under cooperative agreement by Associated Universities, Inc. 
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(2357-114) was observed for 5 minutes after every 30 minutes of time on source. The total 
on-source integration time for WLM in B array was 8.5 hours. The C array observations 
were taken on 2000 March 31 and consisted of a single pointing centered on WLM's position. 
The flux calibrator (0137+331) was observed for 20 minutes in the middle of the observations 
and the phase calibrator (2321-163) was observed for 5 minutes after every 35 minutes on 
source. The total on-source integration time was 5.1 hours. The correlator setup for both 
the B and C array observations was identical to the D array observations; the final channel 
width was 2.57 km s~^ . 

The observations are summarized in Table [TJ Data from each day and each array 
configuration were calibrated sep arately in AIP S using the standard calibration procedure 



outlined in the AIPS Cookbo ok (jGreiserul2002l ). After calibration, the data were imported 



into Miriad (ISault et al.lll995l ) for further flagging as well as continuum subtraction, imaging. 



mosaicking, and cleaning. 

From these data, we produced three final data cubes: two from the D-array mosaic and 
one from the central pointing of the D array mosaic plus the B and C array observations 
(hereafter, the B+C+D array cube). Table [2] gives the final noise and beam size in all three 
cubes. We use the D array mosaic observations to investigate the large scale structure and 
total extent of WLM. After making an initial image to determine the line-free channels, 
the continuum was subtracted from the initial pointings using the task uvlin. Then the 
nine pointings were put together using the mosaic option in the invert task. Two different 
weightings were used: robust equal to 5 (natural weighting) and robust equal to -5 (uniform 
weighting). The resulting data cubes were then cleaned using the mossdi task to deconvolve 
the source distribution from the antenna sampling pattern. Images were cleaned to the 2a 
noise level and the main emission region was selected. Both the MEM and Clean algorithms 
yielded similar results for these observations. The Clean algorithm results were used since 
the final cube had fewer image artifacts after deconvolution. Channel maps for this cube are 
given in Figure [1] and the total intensity (moment zero) image is shown in Figure [2j 

We used the B+C+D array cube to investigate the small scale HI structure of WLM. 
The observations from the B and C configurations were imaged separately using the invert 
task to make sure there were no problems with the data from individual configurations. 
After fiagging, the uv data sets for each array configuration were combined with uvcat. The 
resulting data set was continuum subtracted (uvlin), imaged (invert), and cleaned to the 2a 
noise limit [dean with Miriad choosing the best cleaning algorithm and the main emission 
region selected). The final map was produced using robust equal to zero because this value 
yielded the best beam profile; other values yielded non-Gaussian beam profiles. Figure [3] 
gives channel maps for this data cube. Figures HI [5l and [6] show the total intensity (moment 
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zero) image, the mean velocity (moment one) image, and the velocity dispersion (moment 
two) image for this data cube. 

The moment images were created using the AIPS task MOMNT. We smoothed the 
cubes to three times the resolution in the RA, Dec, and velocity directions. The smoothed 
cube was blanked below the 3a noise limit and used as a mask for unsmoothed cube. The 
channels of the masked, unsmoothed cube are then combined to create a moment zero, one, 
or two image of the data. 

We did not use the combination of the D array mosaic with the central B and C array 
observations in our analysis because the beam (and the noise) would vary significantly across 
the image greatly complicating analysis. A test image created in Miriad combining the D 
array mosaic with the B and C observations does not reveal any new galaxy features that 
are not apparent in the two cubes described above. 



2.2. Optical Images 



Ha images of WLM were obtained using the Perkins 1.8 m telescope at Lowell Obser- 
vatory in 1992 November. The observations used an 800x800 TI CCD provided to Lowell 
Observatory by the National Science Foundation, the Ohio State University Fabry- Perot that 
was used as a 3:1 focal reducer, an Ha filter with a FWHM of 32 A, and an off-band filter 
centered at 6440 A with a FWHM of 95 A. The off-band filter was used to subtract stellar 
continuum from the Ha filter to leave only the Ha nebular emission. In order to survey the 
entire galaxy, two separate pointings were mosaicked to produce the final Ha image. The 
exposure times were 30 00 s. The pixel scale was 0.4 9", and the seeing was ~2.5". Additional 
details can be found in lHunter fc Elmegreenl (120041 ). 



UBV images of WLM were obtained using a SITe 2048x2048 CCD on the Hall 1.1 m 
telescope at Lowell Observatory in three observing runs during 1998 October, 1999 October, 
and 2002 November. The exposures were 3 x 1200 s in V, 4 x 2400 s in B, and 6 x 2400 s 
in U. The seeing was about 3.0", and the pixel scale was 1.13". The telescope position was 
offset between each image in order to allow removal of cosmic rays. The bias was subtracted 
using th e over scan strip, and the images were fiat-fielded using observations of twilight sky. 
Landoltl (119921) standard stars were used to calibrate the photometry. Additional details can 



be found in 



Hunter fc Elmegreenl (120061 ). 
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HI Mass and Extent 



We obtained a total HI flux for WLM of 294 ± 15 Jy km s ^ by integrating the HI spec- 
trum from the naturally weighted D array mosaic cube. See Figure [7] for the spectr um. This 



HI flu x corresponds to an HI mass of 6.3 ± 0.3 x 10^ Mq at distance of 0. 950 Mpc (IDolphin 



2000[) Our flux agrees well with the Parkes Multibeam observations of iBarnes Sz de Blokl 
( 20041 ). which give a total HI flux of 292+^^ Jy km s"^ . Therefore, our interferometric ob 



servations have recovered most of the HI flux in WLM. iJackson et al.l (120041 ) obtained an 
integrated flux density for WLM of only 149 Jy km s~^ from their Australia Telescope Com- 
pact Array (ATCA) observations, which lack the short spacings needed to recover extended 



emission. 



To compare our n aturally weighted, D-arr ay in tegrated intensity map ( Figure [2]) with 
the single dish maps of Barnes &: de Blok (j2004 ) and Huchtmeier et al. J 1981 ), we smoothed 
the D array cube to th e same spatial resolution as Barnes fc de Bloy ( 2004 . 14.4') and 
Huchtmeier et al.l (jl98ll . 9') and created mo ment maps from the smo othe d cubes. Our 



smoot hed maps reach similar noise levels to both lBarnes fc de Blokl (12004 ) andlHuchtmeier et al, 
( 1981 ). We see that our map agrees well with that of Barnes fc de Blok ( 2004h. We do not^ 



howev er, see the irregular , low surface brightness features observed by iHuchtmeier et al. 



( 1l98ll ). IHuchtmeier et al.l (Il98ll ) note that they were not able to separate emission from 
WLM from the nearby Magellanic Stream emission because of their low velocity resolution 
(11 km s~^ ); the tail of th e Magellanic Streani i s very close to the positio n of WLM at an LSR 
velocity of -200 km s~^ ( iPutman et al.ll2003l ). IHuchtmeier et al.l (1 198 if) estimated t he mass 



of WLM by assuming symmetry in WLM's line profile. Although Barnes &: de Blok (j2004 )'s 
spatial resolution is larger ( 14.4' vs. 9'), their superb velocity resolution (1.12 km s~^ , a fac- 
tor of 10 lower than that of IHuchtmeier et al.l ) allows one to easily separate WLM's emission 
from that of the Magellanic Stream in most channels. We, fortun ately, resolve out the Mag - 



ellanic Stream. Given the line profiles shown in Figure [7| and in iBarnes fc de Blo^ 



y (12004) 



(1198 ih 's 



it is clear that WLM's line profile is not symmetric and that IHuchtmeier et al. 
integrated HI mass and maps for this galaxy are not accurate. 

WLM's HI only exten ds to 1.2 times its Holmberg diameter, not the 3.6 times found 
bv IHuchtmeier et al.l (198ll ). We find that the HI diameters quoted for other galaxies in the 
Huchtmeier et al.l (Il98ll ) sample are also t oo large. Deep neutr al hydrogen observations of 



Sextans A with the Green Bank 



gas reported bv IHuchtmeier et al. 



elesco pe (IBowsher et al.ll2003l ) do not detect the extended 
(Il981l) confirming the results of VLA D-array observations 



reported by IWilcots fc Hunten (120021 ). VLA D-array mosaic observations of NGC 2366 re 



duced by Wilcots also do not show signs of the extended HI reported by IHuchtmeier et al. 
(|l98lh . 
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Figure[8]gives the surface density profile for WLM. The inchnation and position angle of 
the galaxy were taken from the natura lly weighted D array r o tation curve (see Section [5]) and 
the disk scale length (i?£)=2.0') from iHunter &: Elmegreen Comparing the surface 

density profile of WLM to the surface density profiles of the other dwarf galaxies in Figure El 
we see WLM's HI distribution is typical of a dwarf irregular galaxy. Several other galaxies 
in this Figure (NGC 2366, NGC 4449, IC 1613, Sextans A, and DDO 88) have central HI 
depressions like WLM's. From Figure [H one can see that the average surface density of 
WLM's ring is about 1.3 times that of the center. The contrast increases to 1.4 times that 
of the center in the portion of the ring with the highest surface density. This ratio is smaller 
than the ratios between ring surface density and central surface density seen in galaxies like 
Sextans A (~ 1.8) and DDO 88 (~ 2) and is more like the ratio seen in DDO 105 (~ 1.2). 
WLM's peak surface density (7.5 Mq pc~^) is average for galaxies in this sample. Its HI 
extends to 8.5 times the optical disk scale length; six of the ten galaxies in Figure [H] have 
larger HI extents. 



4. Small Scale HI Structure 

From Figure HI we can see that the neutral hydrogen distribution of WLM is smooth 
with a larg e central depre s sion a nd slightly asymmetric. Comparing our Figures [2] and H] to 



Figure 4 in IJackson et al.l (120041 ) illustrates why it is crucial to have good uv plane coverage 



for a detailed examination of the HI structure of a galaxy. What we were able to image as a 



broken ring of high column density HI surrounding the center of the galaxy, iJackson et al. 



( I2OO4J ) was only able to image as two asymmetric high column density peaks. We refer to 
the broken ring of high column density HI as the "hook": the handle of the hook is the 
high column density HI extending northward and the hook itself extends from the handle 
south and east around the ring until it ends at the southernmost bit of high column density 
gas. Approximately 20% of the total mass of WLM is in the hook. From the first moment 
image (Figure E]), which shows the velocity centroid distribution, and the position- velocity 
diagram (Figure USD, one can see that the velocity distribution is asymmetric (see Section [S] 
for details). The second moment image in Figure [6] is also asymmetric, with the northern 
half of the galaxy having a higher average velocity dispersion than the southern half. Note 
that the maximum velocity dispersion is approximately 10 km s"^ and the average velocity 
dispersion for the entire galaxy is 4.5 km s~^ . We do not see any sign of the hook expanding; 
plots of spectra in the region as a function of position are single-peaked. 



How does WLM's neutral hydrogen distribution compare to its stellar and ionized hydro- 
gen distributions? Figure [9] shows a color image of the optical data described in Section 12.21 
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with HI contours overlaid. The distribution of the blue stars is closely correlated with the 
high column density HI with the exception of the clump of blue stars on the western edge of 
the galaxy. The HI hook is also well correlated with the Ha emission (see Figure [TUl) except 
for the C2 complex in the southwestern portion of the galaxy. This bubble-like Ha region 
corresponds to an HI deficit at the center of the bubble. Finally, there is faint, diffuse 
Ha emission along the inner eastern edge of the hook indicating that there is on-going star 
formation in this region. 

HI deficits in the center of dwarf irregulars are fairly common (see Figure [HI). These 
rings are thought to be formed by the stellar wi nds and supernova explosions of high mass 
stars (iMcCray fc Kafatoslll987l : lElmegreerull994l ). The high mass stars clear out a bubble in 
the interstellar medium which continues to expand until its expansion velocity slows to the 
velocity dispersion of the galaxy. Gas piles up on the edge of the bubble and can, in turn, 
form more stars. There i s, however , little definitive observational evidence that this process 
is what creates HI rings. lOtt et al.l (120011 ) find evidence that the central ring in Holmberg I 
was the result of a blowout. They find an age for the HI ring that matches the age of the 
stellar populations inside it. Associated hydrodynamic simulations by lVorobyov et al.l (l2004l ) 
show that multi ple supernovae c an pr oduce a feature like the ring seen in Ho I. HST/ACS 
observations by iMomany et al.l (120051 ) of the stellar populations of a portion of SagDIG, 
however, indicate that there were not enough stars formed in th e central region o f the g alaxy 
to generate the observed HI ring. In their study of DDO 88, ISimpson et al.l ( l2005al ) find 
that while the star cluster in the center of the galaxy had enough mass to form the HI ring 
seen today, the cluster is too old and that any HI ring formed by a population as old as the 
cluster would have dissipated by the present day. Is the hole in the center of WLM the result 
of multiple supernova explosions? We can use our observations o f WLM's neutral hydrogen 
along with the star formation properties of WLM (lDolphinll2000l ) to test this hypothesis. 



The first quantity we need to estimate is the timescale for the formation of the hook. 
Assuming that the hook resulted from an expanding shell, we can place a rough upper 
limit on its age. The average central velocity dispersion of WLM is approximately the final 
expansion velocity of the ring if the ring is stalled. Using this velocity (7 km s~^ ) and the 
radius of the hook {Rinner = 250 pc; Router = 462 pc), we obtain an upper limit on the 
expansion time of the ring between 35 Myr (using the inner radius of the ring) and 69 Myr 
(using the outer radius of the ring). An upper limit on the dispersion time for the ring can 
be estimated by finding the time it would take for the ring to expand to twice its observed 
width. Using the same velocity dispersion as above, this time scale is approximately 59 Myr. 
Therefore, the upper limit on the age for WLM's hook (expansion time plus dispersion time) 
is between 94 and 128 Myr. 
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We can use models of a supernova explosion in a uniform medium (jChevalierl 1 19741 ) to 
estimate the energy needed to create WLM's HI hole. To do this, we must first estimate 
the number density of HI at mid-plane (rio). Using an inclination of 75° for WLM and a 
scale height typical of a dwarf irregular galaxy (500 pc), we obtain a pathlength through the 
galaxy of 3800 pc. Dividing the central surface density (1.96 x 10^^ cm~^) by this pathlength 
gives an estimate of no of 0.16 cm~^. A reasonab le upper li r ait on no is 1 cm~^ given the 
properties of HI in other galaxies. Equation (26) of IChevalierl (119741 ) gives an estimate of the 
supernova energy needed to produce a hole of radius R from an explosion into a medium of 
uniform density no. For WLM's HI hole, the total HI kinetic energy is between 3.1 x 10^^ erg 
and 2.45 x 10^^ erg, which corresponds to between 3 and 24 supernovae. 

Has WLM formed enough stars during the last 94 to 128 Myr to produce t hat many su- 
pernovae? Although there have been several stellar population stu dies of WLM (jPerraro et al. 
19891 : iMinniti fc Ziilstralll997l : iReikuba et aPboOol : lDolphin|[200o[ ). there are no spatially re- 
solved star formatio n histories for th e entire galaxy. The highest resolution stellar population 
study was done by iDolphinI (120001 ) . who used WFPC2 on Hubble Space Telescope (HST) 
to get a detailed star formation history for three regions in the western half of the galaxy. 
Figure [TT] shows the location of his fields!^ Unfortunately, his observations missed the high 
column density HI regions, both large HII regions in the galaxy (iHodge fc Millerl (jl995l )'s CI 
and C2 complexes), and the region of blue stars uncorrelated with the hook. In his relatively 
quiescent survey region, he finds that 85% of WLM's star formation in the past 200 Myr is 
in the region surveyed that was closest to the center of the galaxy and that the galaxy has 
been experiencing a burst of star formation in the past Gyr. 

We can see if WLM has had enough recent star formation to form its HI hole. Using a 
typical star formation rate for WLM (1.0 x 10"^ Mq yr~^) scaled to the area interior to the 
rme (1.6 X 10-"^ M0 yr-i) and the maximum lifetime time for the bubble, we derive a total 
mass of s tars forrned be tween 1.5 x 10^ Mq and 2.1 x 10^ M©. Using a Salpeter initial mass 
function (jSalpeter 1955), we find that during the lifetime of the bubble between 138 and 188 
stars were formed with masses greater than 7 Mq, which exceeds the estimated number of 
supernovae needed to form the hole. Therefore, it is plausible that the HI hole in WLM was 
caused by star formation in the center of the galaxy. 

Unlike DDO 88's unbroken ring, WLM's HI ring has a significant gap on the western 
edge of the galaxy and a high column density region extending north along the major axis of 
the disk. The handle of the hook is probably due to an increase in column density because 



^Tlic caption of Figure 1 in Dolphin ( 200d ) is incorrect. The WF2 chip is on the upper left, WF3 is 
the lower left, and WF4 is on the lower right. 
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we are looking through the disk of WLM. The origin of the gap in WLM's HI ring is more 
of a mystery. There are bright, blue stars seen near the C2 Ha complex at the tip of the HI 
hook, which indicates that in the past there was more neutral hydrogen in this region than 
is currently seen. It is likely that the star formation in the C2 Ha region has either used 
up, blown out, or ionized all of the neutral hydrogen. Examining spectra as a function of 
position in the C2 region at the tip of the HI hook, we see double peaked velocity profiles 
characteristic of expanding bubbles with the central channel of the profile being lower than 
the two neighboring channels (see Figure [T^ . The separation between the two peaks is 
approximately 9 km s~^ , which means that the expansion velocity of the bubble is about 
4.5 km s~^ . This velocity is very similar to the velocity dispersion of the galaxy as a whole, 
suggesting that the bubble is ne ar the end of its ex pansion. The Ha emission in this region 
also exhibits signs of expansion (ITomita et al.lll998l . their Figure 5, slit position WLM E). In 
contrast to the C2 region, the current round of star formation in the HI hook (indicated by 
the faint Ha emission along the inner edge of the hook), however, seems to be just beginning 
and it has yet to deplete the neutral hydrogen in this region. The higher velocity dispersion 
in the northern half of the galaxy reflects the disruption of the ISM by the star formation 
taking place in the northern half of the galaxy. A detailed, spatially resolved star formation 
history of the entire galaxy is needed before we can attempt to conclusively correlate neutral 
hydrogen features with star formation episodes. 



Minniti &: Zijlstral (119971 ) have suggested, based on a paper by lWallin et al.l (Il996l ). that 



WLM's globular cluster may have triggered the recent star formation episode in this galaxy a 
Gyr ago. Although we believe that this possibility is unlikely, the globular cluster in WLM, 
assuming a typical globular cluster mass of 10^ M© and a typical WLM halo velocity of 
about 40 km s~^ , does have a total kinetic energy of 1.6 x lO^^erg. If a tenth of this energy 
was transferred to WLM, this energy input would be similar to the energy input by one 
supernova. Without full orbital parameters for the globular cluster, we cannot exclude the 
possibility that a particularly close/violent encounter between WLM and its globular cluster 
has triggered the latest round of star formation. 

Another scenario is that a close encounter with the nearby Cetus dwarf spheroidal 
galaxy may have triggered the star formation. Unfortunately, in this case, we are again 
without detailed orbital parameters. We can constrain the time since interaction, however, 
by assuming that WLM and Cetus dSph are on radial orbits. Then the time since interaction 
is just the distance between the t wo galaxies (200 kpc ) divided by the velocity dispersion of 



the Local Group (~ 60 km s ^ ; Ivan den Berghlll999l ). which yields a time of 3.4 Gyr. This 



time is sar ne order of ma gnitude as the start of the latest burst of star formation in WLM 
a Gyr ago JPolphinlboOol ). 
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Bars can also create a central HI deficit (e.g. iRyder et al.lll996l ). However, we see no 
evidence for a bar in a Spitzer IRAC 3.6 /im image of the central portion of WLM. While the 
approaching side of WLM is clearly warped, we do not see a ny s-shaped isovelo city contours 



in the center of the galaxy that are indicative of a bar (e.g. iRyder et al.lll996l . NGC 1433) 



Rotation Curve 



The position-velocity plot of the B+C+D array data cube in Figure [13] provides a 
first look at the rotation curve of WLM. This plot shows the velocity distribution as a 
function of position along a slice corresponding to the major axis of WLM's neutral hydrogen 
distribution. There are several things to notice about this graph. First, WLM is clearly 
rotating with a velocity of approximately 30 km s~^ . Second, the rotation curve is not 
symmetric. The rotation curve for the approaching side of the galaxy increases to a velocity 
of 30 km s~^ then levels off; the rotation curve for the receding side of the galaxy rises much 
more slowly, does not level off and is slightly distorted at the location of the ring. Finally, 
the ring seen in Figure [His clearly visible in this plot as the two dark spots at ±1.5'. 

To derive a rotation curve for the entire galaxy, we fit a tilted ring model to the velocity 
field of the B+C+D array cube. Tilted ring models divide the galaxy into a series of annuli 
and deterr nine the position angle, inclination, and rotation velocity independently for each 
ring. See iBegemanl ( 19871 ) for more detailed information about this type of model. The 



small beam size of the B+C+D array cube allows us to finely sample the rotation curve and 
minimize beam smearing effects. Initial guesses for the fit parameters were taken from the 
literature. We refined these guesses by fitting functions to the velocity field. To determine 
the center of the galaxy and an approximate inclination, we fit a Brandt type rotation curve 
to WLM's velocity field using the GAL task in AIPS. Then we fit a solid body rotation 
curve to the velocity field to determine the systematic velocity and get a good estimate 
of the position angle. Note that we could not use just one function to get all five initial 
parameters because although the velocity field is better fit by a solid body rotation curve, a 
solid body rotation cu rve cannot determine the center of a galaxy. We used the rotcu r task 
teegeman|[l987l . ll989h in GIPSY JVogelaar fc Terlouwl[200ll : Ivan der Hulst et al.lll992h to fit 
a tilted ring model to the velocity field of the galaxy. The width of each annulus was 15". 
The widths of the annuli were chosen to be slightly less than the major axis of our elliptical 
beam. We started the fits one beam width from the center to avoid beam-smearing effects. 
We have corrected the errors produced by rotcur to take into account that all the pixels in 
the ring are not independent. 



Note that we also fit tilted ring models to the naturally weighted, D array cube velocity 
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fields using the same method. However, since these fits are affected by beam smearing in 
the inner regions and do not yield new information about the rotation curve of WLM, we 
do not discuss them here. 

For the B+C+D array cube, attempting to fit the inchnation, position angle, and ro- 
tation velocity with the same parameters for the entire galaxy led to large residuals with 
spatial structure. Since the approaching (north) and receding (south) halves of WLM are not 
symmetric, we decided to fit the inclination, position angle, and rotation velocity for each 
side of the galaxy separately. The resulting fits are shown in Figure [TH The blue and red 
lines show the fit for approaching and receding sides of the B+C+D array cube, respectively. 
The weighted average of the fitted inclination and position angle for each model are given 
by thin solid lines running across the lower two plots. 

Examining the top panel of Figure [TU we see that the rotation curve for the approaching 
half of the B+C+D array data levels off at about 2.5' to a value of 30 km s~^ , while the 
rotation curve of the receding half continues to rise. The B+C+D array rotation curves 
are relatively bumpy because the B+C+D array cube is sampling velocities on size scales of 
approximately 50 pc, where small scale velocity features from SNe, for example, dominate. 
(See Section m for a discussion of the velocity dispersion distribution of WLM.) 

From Figure [HI one can see that the inclination of the galaxy based on the high reso- 
lution fits is more or less constant as a function of radius. The weighted average inclination 
for each fit is given in Table [31 There is an 8° offset in the inclination between the two sides 
of the galaxy for the B+C+D array rotation curve. 

The position angle as a function of radius shows greater variation than the inclination 
as a function of radius. The position angles of the B+C+D fits change rather abruptly at 
particular radii. In the fit to the approaching side of the B+C+D array cube, the position 
angle changes from about 172° to about 178° at 2.5', decreases to about 170° at 6.5' and 
then begins to increase again at 8'. Note that the position angle for radii greater than 9' is 
very uncertain. The fit to the receding side of the B+C+D array cube increases from about 
170° to about 180° between a radius of 5' and a radius of 8.5'. The average position angle 
for the approaching side is 3° higher than the average position angle for the receding side 
for the high resolution data cube. 

The asymmetry in WLM's rotation curve is not unusual; half the dwarf galaxies in the 
Westerbork HI Surve y of Spiral and Irregular Galaxies (WHISP) are kinemat ically lopsided 



( Swaters et al.l |2002| ) . Only 16% of these galaxies are clearly interacting. ISwaters et al. 



J2002[ ) attributes the rest of the distortions to long-lived remnants of past interactions 



WLM has at least two neighbors that are close enough to have disrupted its neutral hydrogen 
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recently: its ancient globular cluster and the Cetus dwarf spheroidal (~ 200 kpc away). See 
Section H] for a discussion of the interaction of WLM with these objects. The presence of 
a globular cluster suggests that WLM interacted with something a long time ago; many of 
the dwarfs that have young s uper star clusters in them today have nearby HI clouds or show 
other signs of an interaction (IBillett et al.ll2002l ). 

We can use our derived rotation curves to fit a mass model to the data and determine 
the dynamical mass of WLM. The details of this procedure are in Appendix |X1 We derive a 
dynamical mass of 2.16 x 10^ Mq for WLM. Our mass models show that WLM is ex t remel y 
dark matter dominated when compared to the dwarf irregular galaxies in ICote et al.l (120001 ) . 



Jackson et al.l (120041 ) derive d a rotation curve for WLM from their low resolution HI 



data. Their rotation curve (see iJackson et al.l |2004| . Figure 8) is very similar to our curve 
(see our Figure [T^ . It reaches a velocity of about 38 km s~^ at a radius of 4 kpc, while 
our low resolution curve reaches a velocity of about 42 km s~^ at the same radius. We also 



derive similar halo dynamical massesjfl 



6. Summary and Conclusions 

Dwarf irregular galaxies are unique laboratories for studying the interaction between 
stars and the interstellar medium in low mass environments. In this paper, we have presented 
the highest spatial resolution observations to date of the neutral hydrogen content of the 
Local Group dwarf irregular galaxy WLM. These observations have allowed us to study the 
interstellar medium of WLM in detail. The major conclusions of this work are: 



• WLM's neutral hydrogen distribution is typical for a galaxy of its type and size. Our 
synthesis observations were able to recover all of the HI flux from this galaxy; we derive 
an HI mass of 6.3 ± 0.3 x 10^ M^:, for WLM. We derive an HI extent of ~ 30', which 



is much less than the 45' extent derived by iHuchtmeier et al.l (Il98ll ). In the central 
region of WLM, the HI has a hook- like distribution, i.e., a central broken ring of HI 
with a linear feature extending northwards from the top of the broken ring. 

• We shoe that the hook is likely the result of star formation propagating out from 
the center of the galaxy. The young stars and Ha emission in this galaxy are mostly 
spatially correlated with the neutral hydrogen peaks. The gap on the western edge of 



■^Note that the dynamical mass quoted in iJackson et al.l (j2004l ) is incorrect. According to D. Jackson 
(2006, private communication), the correct dynamical mass is 1.4 x 10^ M© . 
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the hook is likely the result of star formation in that region using up, blowing out, or 
ionizing all of the neutral hydrogen. The diffuse Ha on the eastern edge of the hook 
indicates that there is ongoing star formation in the central region of WLM. 

• WLM's velocity field is asymmetric. The approaching (northern half) appears to be 
warped and the inner approaching region has a steeper velocity gradient than the inner 
receding region. We derive a dynamical mass for WLM of 2.16 x 10^ Mq. 

In short, WLM is a typical dwarf irregular galaxy undergoing a recent burst of propa- 
gating star formation that is disrupting its neutral hydrogen. A detailed, spatially resolved 
star formation history of the entire galaxy would allow us to track the progress of the star 
formation across the galaxy and attempt to correlate particular HI features with specific 
episodes of star formation. 
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A. Deriving the Dynamical Mass of WLM 

We can use our derived rotation curves to fit a mass model to the data and de- 



associated with fitting mass models (see 


Swaters et al. 


( 

); 


2000h: 


de Blok et al. 


(2001a 


b); 


van den Bosch & Swaters ( 


2001) 


: de Blok & Bosma 


c 


2002 


de Blok et al. 


(2003 


); 


VIcGaueh et al. 


(2003 


); 


Swaters et al. 


f2003aljbh: 


Weldrake et al. 


(200: 


5); 


de Blok 


(2005 


); 


Gentile et al. 


(2005) 



for a sampling). Here we outline our procedure. 

For each side of WLM, we used a combination of the B+C+D array and the D array 
mosaic rotation curves to fit a mass model. The D array mosaic rotation curves were gen- 
erated the same way as the B+C+D array rotation curves (see Section Ej). The D array 
mosaic rotation curve extends to larger radii than the B+C+D array rotation curve, which 
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is crucial for constraining the halo model. The B+C+D array rotation curve is less affected 
by beam smearing in the inner regions than the D array mosaic rotation curve, so it is better 
at constraining the disk model. We spliced the B+C+D rotation curve to the D array mosaic 
curve at 9' for the approaching side and at 5' for the receding side. The B+C+D array and 
the D array mosaic rotation curves agree in the overlap region. 

We used a simple model consisting of a gas disk, a stellar disk, and a halo to determine 
the dynamical mass of WLM. For the gas disk, we found the neutral hydrogen contribution 
to the total rotation curve from the HI surface density profile, whose inclination and position 



Beeeman 


(1987 


) and 


Casertano ( 


1983) 



(which is proportional to the HI mass distribution) was then multiplied by 1.43 to account 
for the presence of helium and other gas. The contribution of the opti cal component to the 



total rotation curve was found from a V-band surface brightness profile (IHunter fc Elmegreen 



20061 ) using the same method. We converted the luminosity distribution of the galaxy into a 
mass distribution (which is proportional to the rotation velocity squared) by multiplying the 
luminosity distribution by a mass-to-light ratio. The mass-to-light ratio is a free parameter 
in our fits. It was restricted to physically meaningful values, i.e., the mass-to-light ratio had 
to be greater than zero. We also assumed that the mass-to-light ratio did not change as a 
function of radius in the galaxy. For the halo , we used two different models: an isothermal 
sphere and a NFW halo ( INavarro et al.lll997l ). The first model is generally a good fit to the 
data and the second model is suggested by cosmological simulations. The circular velocity 
for the isothermal sphere with density distribution 



is given by 



p(r) = po [l + ir/af] 



AnGpoo^ll — (a/r) arctan(r/a)]. 



(Al) 



(A2) 



wher e p is the density, r is the r adius, po is the central halo density, and a is the core radius 



[e.g. 



McGaugh &: de Bloklll998l ). For the NFW halo, the circular velocity is given by 



^200 



ln(l + cx) — cx/(l + cx) 



-^c,nfw 



ln(l + c) -c/(l + c) 



X 



(A3) 



where r2oo is the radius where the overdensity is 200 times that of the background density, 
f2oo is the velocity at r2005 ^ is 'r/r2Qo, and c is the concentration parameter. The parameters 
r2oo and f 200 are related by the equation 



V200 



GM 
^200 



1/2 



^200 
h^^ kpc 



1/2 



:i + zof^ km s-^ 



(A4) 



- 16 - 



where Qq is the current matter density and ^(2:0) is the matter density at redshift zo- The 
redshift of WLM was taken to be 0. The concentration parameter is just r2oo divided by the 
scale radius of the halo so it was restricted to be greater than or equal to one (r2oo > ''^s)- 

We fit three different types of mass models to the data: a maximum disk model, a 
maximum halo model, and a best fit mass model. The maximum disk model fits the data 
assuming a maximal contribution from the disk to the rotation curve. Note that this mass 
model still includes a halo component to the fit since a disk rotation curve cannot replicate 
the fiat portion of the observed rotation curve. The maximum halo model fits the data 
assuming that there is no contribution to the rotation curve from the stellar or gas disk. 
These two fits, while clearly not physical, provide important boundary conditions for our 
mass model of WLM. The best fit mass model is a model of the galaxy which includes 
contributions from the halo, the stellar disk, and the gas disk. For this fit, the parameters 
for each component of the mass model are chosen to minimize the difference between the 
observed and the model rotation curves using the parameter as our m inimization criteri on. 



We used the downhill simplex (amoeba) method of Nelder and Mead (jPress et al.lll992l ) to 
find the mass model parameters which minimized the value. While this might not be the 
best algorithm for minimizing the parameter, it does have the advantage of being easy 
to implement and relatively fast given that we only had to fit three parameters. To ensure 
that the algorithm found a truly global minimum, we ran it with a grid of different initial 
guesses for the parameters and chose the solution that had the lowest chi-squared value. See 
Tables S] and E] for a summary of the mass model fits and Figure [TS] for plots of the best fit 
mass models. 

We were not able to generate an acceptable mass model fit to the entire rotation curve 
for the approaching side. The rotation curve for the approaching side (see the plots on the 
left hand side of Figure [T5I) increases rapidly at radii greater than 8', deviating from the 
behavior of a well-behaved rotation curve (cf. the rotation curve for the receding side of the 
galaxy in show in the plots on the right hand side of Figure [T5l) . The warp can also be seen 
in Figure [51 The isovelocity contours form a U-shape on the northern (approaching) side of 
the galaxy, while the contours on the southern (receding) side remain more or less straight. 
Unfortunately, there is not enough information in either the B+C+D rotation curve or the 
D array mosaic rotation curve to constrain the change in inclination and/or position angle. 
However, the unusual rotation curve generated for the approaching side obtained by fixing 
both the inclination and position angle to a single set of values points to unconstrained 
variations in the inclination and position angle. Since the approaching side of the galaxy is 
affected by a warp at radii greater than 8', we excluded these radii from our fit. Note that 
the results in Tables H] and [5] for the approaching side only go out to the last point we fitted. 
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Due to the apparent warp in the approaching side of WLM, we used the isothermal 
sphere best fit mass model to the receding side of the galaxy to determine an estimate for 
the dynamical mass of the galaxy (2.16 x 10^ M©) and the mass model parameters for the 
entire galaxy. The maximum halo model (no disk) for the receding side yields a relatively 
good fit to the data for the isothermal sphere case. The maximum disk model, however, yields 
poor Xr values for all fits. WLM is clearly not a galaxy whose dynamics are dominated by its 
disk. Note that WLM's estimated total dynamical mass changes by less than 10% between 
the isothermal sphere halo and the NFW halo best fits for the undisturbed receding side of 
the galaxy. 

We also compared our best fit mass model (isothermal sphere halo plus gas and stellar 
dis ks) parameters for WLM with those derived for a sample of five dwarf irregular galaxies 
by ICote et al.l (120001 ). We find that the core radius and central density for the receding 
side of the rotation curve as well as the derived galaxy parameters (stellar M/ L ratio, total 
dynamical mass, and ratio of mass in dark matter to the mass in luminous matter) are within 
the range of values established b y these auth ors. The halo parameters for the approaching 
side of the galaxy fall outside the ICote et al.l distribution, but these parameters are not well 
constrained for th e the appro aching side because of the warp. WLM, when compared with 
the galaxies in the ICote et al.l sample, is on the extremely dark matter-dominated end of the 
dwarf galaxy spectrum. 
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Table 1. Summary of Observations 



Array Configuration B CD 

Date 2002 July 9 and 12 2000 March 31 2001 October 19 

Number of Pointings 1 19 

Square Field of View (') 30 30 60 

Channel Width (km s"^ ) 2.57 2.57 2.57 

Flux Calibrator 0137+331 0137+331 0137+331 

Phase Calibrator 2357-114 2321-163 2321-163 

Total Integration Time (hours) 8.5 5.1 ~2 



Table 2. Summary of Final Data Cubes 









Beam Size 


1(7 Noise per Channel 




Cube Name 


Weighting 


(") 


(mJy beam~^) (K) 


D 


array mosaic 


Natural 


83.84 X 59.11 


4.7 0.58 


D 


array mosaic 


Uniform 


72.77 X 48.88 


5.1 0.87 


B- 


f C+D array cube 


Robust=0 


16.17 X 10.59 


1.4 5.0 



Table 3. Tilted Ring Rotation Curve Fit Results 



Data Set 


Fit 


<P.A.> 

n 


n 


<Incl> 

n 


n 


D array mosaic 


Entire Galaxy 


179.3 


2.2 


69.8 


3.7 


D array mosaic 


Approaching 


174.6 


2.3 


62.8 


6.0 


D array mosaic 


Receding 


179.5 


2.0 


72.2 


3.2 


B+C+D array 


Approaching 


175.5 


3.8 


74.2 


3.1 


B+C+D array 


Receding 


172.3 


4.1 


82.4 


5.8 
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Table 4. Summary of Mass Model Fits with Isothermal Sphere Halos 





a 


Po 


M/L 


Mgas 


-^stars 


Mhalo 


Fit Type 


n 


(Mo pc-3) 


(Mo /Lo) X'r 


(Mo) 


(Mo) 


(Mo) 




A]:)]:)ro aching Side 


Best Fit 


71.5 


0.157 


2.00 


6.46 X 10^ 





3.82 X 10^ 


Max. Halo 


83.3 


0.150 


2.57 








4.75 X 10^ 


Max. Disk 


0.00985 


174000 


4.8 9.90 


6.46 X 10^ 


2.07 X 10^ 


1.02 X 10^ 


Receding Side 


Best Fit 


337 


0.0265 


1.1 0.98 


9.91 X 10^ 


4.73 X 10^ 


2.01 X 10^ 


Max. Halo 


237 


0.0483 


0.84 








2.15 X 10*^ 


Max. Disk 


6870 


0.00737 


3.6 8.88 


9.91 X 10^ 


1.55 X 10^ 


2.61 X lO'^ 



Table 5. Summary of Mass Model Fits with NFW Halos 









M/L Mgas 


■^stars 


Mhalo 


Fit Type 


(km s~^ ) 


c 


(Mo /Lo) xl (Mo ) 


(Mo) 


(Mo) 




Approaching Side 


Best Fit 


24.6 


6.60 


0.965 2.34 6.46 x 10^ 


4.17 X 10^ 


3.24 X 10^ 


Max. Halo 


32.1 


6.81 


5.42 





5.07 X 10^ 


Max. Disk 


3.33 


2.28 


4.8 10.42 6.46 x 10^ 


2.08 X 10^ 


5.82 X 10^ 


Receding Side 


Best Fit 


149.6 


1.0 


7.73 9.91 X 10^ 





2.23 X lO'^ 


Max. Halo 


165.4 


1.0 


8.55 





2.50 X 10*^ 


Max. Disk 


87.3 


1.0 


3.6 63.94 9.91 x 10^ 


1.56 X 10^ 


1.26 X lO'^ 
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Fig. 1. — Emission in each channel for the naturally weighted D array mosaic of WLM. 
Every other channel is shown. The contours are 1, 2, 5, 10, 20, and 30 times the 3cr noise 
level (14.1 mJy beam~^; 0.56K). The beam is given in the lower left hand corner of the 
first panel. The channel maps for the uniformly weighted D array mosaic of WLM are very 
similar. 



Fig. 2. — Integrated intensity map of WLM from the naturally weighted D array mosaic. 
The contour levels are 1, 10, 20, 30, 50, 100, 200, and 300 times the 3a noise level (8.1x10^^ 
cm~^). The greyscale image is provided to help distinguish the contours. The emission at the 
edge of the plot is result of the mosaicking process (we have better sensitivity at the center 
of the image due t o the overlap of pointings than at the edges) and is not real emission. 



Barnes fc de Blokl (120041 ) find no companion HI clouds near WLM with a detection limit of 



1.9 X lO^M© for a object with a velocity dispersion of 10 km s ^ . The beam is given in the 
lower left hand corner of the image. 



Fig. 3. — Emission in each channel for the B+C+D array cube. Every other channel is 
shown. The baryocentric velocity of the channel is in the upper right of each panel. The 
contours are for 3, 10, and 20cr emission. The la noise is 1.4 mJy beam~^, which is equal to 
5. OK for our beam. The beam is shown in the last panel. 



Fig. 3. — Continued 



Fig. 4. — Integrated intensity (moment zero) map of WLM from the B+C+D array cube. 
The contour levels are 3, 20, 45, 90, and 120 times the la noise level (2.9x10^^ cm~^). The 
greyscale image is provided to help distinguish the contours. The beam is shown in the lower 
left hand corner of the image. 



Fig. 5. — Mean velocity map (moment one) of WLM from the B+C+D array cube. The 
contour levels are -150, -140, -130, -120, -110, -100 and -90 km s~^ . The greyscale image is 
provided to help distinguish the contours. The beam is shown in the lower left hand corner 
of the image. 



Fig. 6. — Velocity dispersion (moment 2) map of WLM from the B+C+D array cube. The 
contour levels are two and three times the width of one channel (2.57 km s~^ ). The greyscale 
image is provided to help distinguish contours. The beam is shown in the lower left hand 
corner of the image. 



Fig. 7. — Integrated HI spectrum of WLM. The errors in this spectrum are on the order of 
the width of the plotted line. 



Fig. 8. — Surface density profile for WLM (black dots and black triangles connected with 
a solid black line) compared to the surface density profiles of other dwarf galaxies (col- 
ored lines). The black triangles show the surface density profile for the B+C+D array 
cube, while the black dots show the surface density profile for the naturally weighted D 
array mosaic cube. Rd is the V band disk scale length from iHunter fc Elmegreeru (120041 ). 
For WLM, = 2. 0^ The refer e nces for the surfac e density prof il es of the other galax- 
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and DDO 43, ISimpson et al.l fl2005bf ). 



Fig. 9. — A combined color image of the UBV optical data with the HI contours from 
Figure H] overlaid. 



Fig. 10. — HI contours overlaid on the Ha image of WLM. The contour levels are 1, 30, 50, 
70, 1 00, 130, 160 tim e s the la noise level (2.9 xlO^^cm'^). The CI and C2 Ha complexes 
from iHodge fc Milled (119951 ) are indicated by arrows. 



Fig. 11.— The location of the lDolphinI (120001 ) HST WFPC2 fields overlaid o n a DSS image of 
WLM . The red squares indicate the position of the Ha regions identified by lHodge &: Miller 
(]l995h . The CI and C2 complexes are direc tly north and south o f the easternmost WFPC2 
field. This image was created using Aladin (jBonnarel et al. 2000). 



Fig. 12. — HI spectra for each pixel in the C2 Ha complex. Note the double peaked profiles 
near the center of the region. The separation between the two peaks is approximately 9 
km s~^ , which gives an expansion velocity for the bubble of 4.5 km s~^ . 

Fig. 13. — Position velocity diagram for a cut along the major axis of the neutral hydrogen 
distribution of WLM (PA = 179.6°, which is the position angle derived for the D array 
mosaic in Section [5]). Positive offsets are along the receding axis and negative offsets are 
along the approaching axis. 
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Fig. 14. — Summary of tilted ring model rotation curve fits. The top panel shows the 
rotation velocity as a function of radius, the middle panel shows the inclination as a function 
of radius, and the bottom panel shows the position angle as a function of radius. The two 
fits are shown in different colors: blue lines are for the fit to the approaching side of WLM 
using the B+C+D array cube and red lines are for the fit to the receding side of WLM using 
the B+C+D array cube. Error bars for each curve are shown in the corresponding color. 
The error bars for some points are quite large and extend vertically all the way across the 
graph. The thin solid lines running across the inclination and position angle plots are the 
weighted average of the inclinations and position angles for each fit. 



Fig. 15. — Best fit mass models for WLM. The left hand plots are for the approaching side 
of the galaxy and the right hand plots are for the receding side of the galaxy. The top row 
shows the fit for an isothermal sphere halo and the bottom row shows the fit for an NFW 
halo. The rotation curve is a combination of the B+C+D array and the D array mosaic 
rotation curves. The first two points were excluded from the fits to either side because they 
are the points most affected by inclination and beam smearing effects. Points at radii greater 
than 8' were excluded from the fit to the approaching side because beyond this radius the 
galaxy is warped and we are unable to constrain its inclination. The dashed line is the 
contribution to rotation curve from the best fit halo, the dotted line is the HI contribution, 
the dash-dot curve is the total disk contribution, and the solid line is the total rotation curve 
(disk plus halo). Note that because the derived M/L ratio of the disk is zero for some fits, 
the HI contribution and the total disk contribution lines overlap and there is no plotted 
contribution from the stellar distribution. 
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